AD-750  540 

CODA  SUPPRESSION  CAPABILITIES  OF  THE  BEAM 
AND  MIXED-SIGNAL  PROCESSOR 

T.  J.  Cohen 
Teledyne  Geotech 


Prepared  for: 

Advanced  Research  Projects  Agency 
30  August  1972 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


BEST 

AVAILABLE  COPY 


AD750540 


298 


CODA  SUPPRESSION 
CAPABILITIES  OF  THE  BEAM  AND 
MIXED-SIGNAL  PROCESSOR 


T.J  COHEN 
SEISMIC  DATA  LABORATORY 

30  AUGUST  19/2 


Pvtptiwd  ftr 

AIR  FORCE  TECHNICAL  APPLICATIONS  CENTER 

Washington,  D  C. 

Project  VELA  UNIFORM 

Sa»n»rti  tf 

ADVANCED  RESEARCH  PROJECTS  AGENCY 
NuJeor  Monitoring  Research  Office 
ARPA  Order  No.  1714 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

1  S  L)»[)or'i  «»  »  f  C 

St  r  pH  J  A  J2  SI 

ALEXANDRIA  LABORATORIES 


^TELEDYNE  GEOTECH 


wue-sa  rente  Btu-ntt; 


"STiimioiimiHjrEB. 


Security  Certification  _ —  ** 

DOCUMENT  CONTROL  DATA  •  *10  „  „  .  .... 

„ — r  ......  - - *  ~  .7".".-. : — *sr 

I  o«ioiu»tin  c  actwity  rc.*w«f«  •*<»**)  .  , 

Teledvne  Geotech  Unclassified - 

Alexandria,  Virginia  *'  •"00’, 


I  «»0«T  TlTLl 

CODA  SUPPRESSION  CAPABILITIES  OF  THE  BEAM  AND  MIXED-SIGNAL 

PROCESSOR 

Scientific  _ 


t  AUTMOW*)  (L—i  Urn*  mill 

Cohen,  T.J. 


a.  mtPonr  date 

30  August  1972  _ _____ 


to  CONTRACT  ON  tWAMT  NO- 

F33657-72-C-0009 

fc  roojicT  no 

VELA  T/2706 

ARPA  Order  No.  1714 
rfARPA  Program  Code  No.  2F-10 


,0  AVAIL  ABILITY /LIMITATION  NOTICII 

Wf MW  HI  M»uc  mtm; 

[Advanced0 ReseaVch V*Pro j  ects  Agency 
Nuclear  Monitoring  Research  Office 
Washington,  D.C. _ _ 

Usifci  7.  and  19-eieaent  subarray;  at  TFO  wa  coapara  tha  relativa  ibiltttas  of  the  baaa  and  a 
aixad-signai  procassor  to  attenuate  P-wave  codas.  P-wave  signals  froa  va/ious  earthquakes  ara 
tiae-ihifted  to  simulate  arrivals  froa  diffarant  azlauths  and  6C’  dlstanca.  Thase  signals  are 


ft.  giTM<EW^wyP°WT  MOfi)  f4 1»r  iwwaWw  mMf 


»1  #‘a*T*ACT 


to  axcaed  that  of  tha  beaa.  Up  to  14  ob  laproveaant  uver  tha  baaa  Is  obtainad,  although  taprove- 
aent  is  ganeraily  on  tha  order  of  3  to  S  db  for  both  subarrays.  A  noainal  valua  of  11  db  Is 
repreventat iva  of  tha  aaiiaua  coda  attanuation  obtained  using  tha  Mixed  signal  procassor  and  tie 
19-e laaent  subarray:  for  tha  7-eiaaent  subarray,  aaxlaua  attenuation  is  roughly  14  db.  Coda 
attanuation  obtained  using  the  7-eieaent  subarray  and  tha  alied-signal  processor  is  coaparable  to 
that  obtatned  using  the  19-eiaaant  subarray  and  tha  beaa.  Representative  TFO  coda-attanuatlon 
curves  for  tha  alxad-signal  processor  ara:  19  eleaents: 

S.74X  •  0.91X2 

o  •  - — - y  dh 

i.O  •  0.31X  •  0.043X* 


u— 2  ‘JUL*  0-69X* 

1.0  ♦  0.53X  ♦  O.00SX* 


*  ’  1 ,  '  'al’j1 

For  X  snail,  the  aiied*  s  ignal  processor  is  sis  tlaei  aore  sensitive  to  chan,et  in  X  than  it  the  beta. 
P.ellainary  remits  suiaett  th.t  uhen  the  difference  in  the  ray  parameter  vectors  is  ta.ll,  a 
l.rie-.perture  array  and  the  aiaed*sl|nal  processor  are  required  for  si,nificant  coda  attenuation 
(-9  db). 


!«  KEY  WOROS 


Beamforming 

Mixed-Signal  Processing 
Coda  Suppression 


9 


CODA  SUPPRESSION  CAPABILITIES  OF  THE 
BEAM  AND  MIXED-SIGNAL  PROCESSOR 

SEISMIC  DATA  LABORATORY  REPORT  NO.  298 

VELA  T/2706 


AFTAC  Project  No.: 
Project  Title: 

ARPA  Order  No. : 

ARPA  Program  Code  No.: 

Name  of  Contractor: 

Contract  No. : 

Date  of  Contract: 

Amount  of  Contract: 
Contract  Expiration  Date: 
Project  Manager: 


Seismic  Data  Laboratory 

1714 

2F-10 

TELEDYNE  GEOTECII 

F33657-72-C-0009 

01  July  1971 

$  1,736,000 

31  October  1972 

Robert  R.  Blandford 
(703)  836-3882 


P.  0.  Box  334,  Alexandria,  Virginia  22314 


II 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


TABLE  OF  CONTENTS 


Page  No. 

ABSTRACT 

INTRODUCTION  1 

SIGNAL  ANALYSIS  THEORY  3 

DATA  ANALYSIS 

o 

RESULTS 

CONCLUSIONS  15 

REFERENCES 

ACKNOWLEDGEMENTS  ,  0 


IK 


IMMU 


LIST  OF  FIGURES 


Figure  Title  Figure 

Least-square  technique  for  estimation  of 

one  signal  in  presence  of  another  (example).  1 

Geographical  location  of  TFO.  2 

Array  configuration  for  TFO.  3 

Seismograms  for  the  Easter  Island  event.  4 

Seismograms  for  the  Fiji  islands  event.  5 

Seismograms  for  the  Fox  Islands  event.  6 

Seismograms  for  the  Hokkaido  event.  7 

Seismograms  for  the  Tonga  Islands  event.  8 

Coda  attenuation,  19-element  subarray, 

Easter  Island  event,  seismograms  bandpass 
filtered  with  3  db  points  at  0.4  and  3.C  Hz.  9 

Coda  attenuation,  7-element  subarray,  Easter 

Island  event,  seismograms  bandpass  filtered 

with  3  db  points  at  0„4  and  3.0  Hz.  10 

Coda  attenuation,  19-element  subarray,  Fiji 
Islands  event,  seismograms  bandpass  filtered 
with  3  db  points  at  0.4  and  3.0  Hz.  11 

Coda  attenuation,  7-element  subarray,  Fiji 
Islands  event,  seismograms  bandpass  filtered 
with  3  db  points  at  0.4  and  3.0  Hz.  12 

Coda  attenuation,  19-element  subarray,  Fox 
Islands  event,  seismograms  bandpass  filtered 
with  3  db  points  at  0.4  and  3.0  Hz.  13 

Coda  attenuation,  7-element  subarray,  Fox 

Islands  event,  seismograms  bandpass  filtered 

with  3  db  points  at  0.4  and  3.0  Hz.  14 


LIST  OF  FIGURES  (Cont'd.) 


Figure  Title 

Figure  1 

Coda  attenuation,  19-element  subarray, 
Hokkaido  event,  seismograms  bandpass 
filtered  with  3  db  points  at  0.4  and 

3.0  Hz. 

15 

Coda  attenuation,  7-element  subarray, 

Hokkaido  event,  seismograms  bandpass 
filtered  with  3  db  points  at  0.4  and 

3.0  Hz. 

16 

Coda  attenuation,  19-element  subarray, 

Tonga  Island  event,  seismograms  bandpass 
filtered  with  3  db  points  at  0.4  and  3.0  Hz. 

17 

Coda  attenuation,  7-element  subarray,  Tonga 
Island  event,  seismograms  bandpass  filtered 
with  3  db  points  at  0.4  and  3.0  Hz. 

18 

Coda  attenuation,  19-element  subarray.  Fox 
Islands  event,  seismograms  bandpass  filtered 
with  3  db  points  at  0.4  and  3.0  Hz. 

19 

Coda  attenuation,  7-element  subarray,  Fox 
Islands  event,  seismograms  bandpass  filtered 
with  3  db  points  at  0.4  and  3.0  Hz. 

20 

Coda  attenuation,  19-element  subarray,  Tonga 
Islands  event,  seismograms  bandpass  filtered 
with  3db  points  at  0.8  and  2.0Hz. 

21 

Coda  attenuation,  7-element  subarray,  Tonga 
Islands  event,  seismograms  bandpass  filtered 
with  3db  points  at  0.8  and  2.0Hz. 

22 

Average  coda  attenuation,  19-element  sub¬ 
array  . 

23 

Average  coda  attenuation,  7-element  sub¬ 
array. 

24 

Attenuation  as  a  function  of  the  difference 
in  ray  parameters,  19-element  subarray,  beam 

25 

LIST  OF  FIGURES  (Cont'd.) 


Figure  No. 


Figure  Title 

Attenuation  as  a  function  of  the 
difference  in  ray  parameters, 

19-element  subarray,  mixed-signal 
processor. 

Attenuation  as  a  function  of  the 
difference  in  ray  parameters,  7-element 
subarray,  beam. 

Attenuation  as  a  function  of  the  differ¬ 
ence  in  ray  parameters.  7-element  sub¬ 
array,  mixed-signal  processor. 

Coda  attenuation,  19-element  subarray, 
Fox  Islands  event,  seismograms  bandpass 
filtered  with  3  db  points  at  0.4  and 
3.0  Hz.  (Signal  2  from  Tonga  Islands). 

Coda  attenuation,  7-element  subarray, 

Fox  Islands  event,  seismograms  bandpass 
filtered  with  3  db  points  at  0.4  and 
3.0  Hz.  (Signal  2  from  Tonga  Islands). 

Comparison  of  coda  attenuation  figures, 
Easter  Island  event,  19-element  subarray 
and  wide-aperture  7-element  subarray. 

Comparison  of  coda  attenuation  figures, 
Easter  Island  event,  small-aperture 
7-element  subarray  and  wide-aperture 
7-element  subarray. 


26 


27 


28 


29 


30 


31 


ii 


list  of  tables 


Table  Title 

Earthquake  Source  Data 

Coordinates  of  TFO  Instruments 

Coda  Attenuation  Parameters 

Coda  Attenuation  Parameters 

Coda  Attenuation  Parameters 

Comparison  of  Coda  Attenuation 
Capabilities 


Table  No. 

I 

II 

III 


IV 

V 

VI 


INTRODUCTION 


The  number  of  opportunities  for  masking  the  signals 
from  an  explosion  in  the  coda  of  an  earthquake  depends 
critically  on  the  amplitude- time  characteristics  of  the 
earthquake  coda  at  various  stations  within  a  monitoring 
network.  When  coda  are  suppressed,  the  time  window 
available  for  masking  is  reduced,  and  so,  correspondingly, 
are  the  number  of  testing  opportunities. 

Given  that  the  explosion  is  to  be  masked  in  an 
earthquake  which  may  occur  anywhere  on  the  earth,  the 
classical  method  for  coda  suppression  is  to  beam  on 
the  suspected  test  site.  For  a  19-element  subarray  such 
as  the  innermost  elements  at  TFO ,  bermforming  might  be 
expected  to  yield  coda  suppression  on  the  order  of  13  db. 
This  attenuation  figure,  however,  is  based  on  a  random 
noise  model.  In  reality  the  noise  we  are  trying  to 
suppress  is  a  teleseismic  signal  propagating  across  the 
array.  We  therefore  expect  coda  attenuation  to  be  a 
function  of  the  difference  between  the  apparent  veloci¬ 
ties  across  the  array  for  the  two  interfering  signals, 
as  well  as  of  the  angular  separation  between  the  event 
epicenters.  Further,  because  a  mixed-signal  processor 
would  take  account  of  a  second  signal's  presence,  we 
expect,  a  priori,  that  this  processor  would  equal  or 
exceed  the  coda  attenuation  capability  of  the  beam. 

In  this  report  we  investigate  the  coda  suppres¬ 
sion  capabilities  of  the  beam  and  mixed-signal  pro¬ 
cessor  using  7-  and  19-element  subarrays  at  TFO.  The 
analysis  consists  of  determining  how  much  of  signal  1 
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(the  earthquake)  "leaks"  into  our  estimate  of  signal  2 
(the  event  to  be  masked).  We  would  like  to  perform  this 
analysis  using  two  overlapping  signals  actually  recorded 
at  TFO.  However,  this  would  not  yield  acceptable  results 
because  we  would  not  know  the  true  waveforms  of  either 
signal.  Alternatively,  we  could  add  two  signals  together. 
In  this  case,  we  would  have  to  form  the  beam  of  the 
signal  from  the  explosion  site  and  subtract  it  from  the 
beam  of  the  composite  seismograms  in  order  to  determine 
the  residual  signal  leaking  through  from  the  earth¬ 
quake  signal.  This  method  would  have  to  be  employed  if 
nonlinear  processes  were  involved.  Fortunately,  both 
the  beam  and  mixed-signal  processor  involve  linear 
processes.  Thus,  signal  2  need  not  be  present  in  the 
seismograms  analyzed  because  we  only  seek  to  determine 
how  much  of  signal  1  is  present  in  our  estimate  of 
signal  2.  We  proceed  as  follows:  we  simply  simulate  the 
arrival  of  an  earthquake  at  various  azimuths  and  at  a 
fixed  distance  of  60°  from  TFO,  and  estimate  the  signal 
leakage  for  a  fixed  test  site.  The  ratio  of  the  root 
mean  square  (rms)  power  for  the  leakage  signal  to  the 
average  rms  power  of  the  earthquake  arrivals  is  then 
defined  as  the  coda  attenuation. 

We  analyzed  the  seismograms  for  five  real  events, 
time-shifted  them  to  produce  simulated  arrivals,  and 
studied  the  coda  attenuation  for  a  variety  of  signal 
waveform  types. 


-2- 


SIGNAL  ANALYSIS  THEORY 


Consider  the  following  representation  for  a 
collection  of  signals: 

oo  P 

=  l  l  xik(t-u)s.(u)  +  n .  (t) 

J  u=-°°  k  =  l  J K  K  J 

where 


y j  J  =  !»•••»  ti  is  the  collection  of  N  observed 
series ; 

s(t)»  k  -  1,...,  p  is  the  collection  of  signals  to 
be  estimated 

xjk^t^  =  6^t-Tjk^  w^cre  6(t)  *  1  for  t  -*  0,  zero 
otherwise] ; 

and 

nj (O  is  th^  noise  on  the  j'th  channel. 

Further,  assume  that  the  signals  can  be  estimated  by  a 
linear  estimates  of  the  type: 

oo 

§k(t)  =  .1  J  hkj ( x  ^  y  j (t-x)dx 

J  ®  1  m  OO 

where  h^  (t)  is  a  p  x  n  matrix  of  filter  functions  to 
be  determined. 


Using  the  above,  it  can  be  shown  (Dean  et  al., 

1968)  that  in  the  frequency  domain  the  matrix  of  filter 
coefficients  H(w)  which  produces  the  best  linear  unbiased 
estimates  of  the  sR(t)  is  given  by 

u(u>)  = 


where  an  n  x  n  matrix  of  spectra  and  cross¬ 

spectra  at  frequency  u.  By  calculating  11(e))  at  each 
frequency  and  then  transforming  back  into  time,  one 
need  only  convolve  the  filters  obtained  with  the  y.{t) 
to  produce  the  signal  estimates  sk(t).  lor  a  single 
signal  (k=l)  in  unrorrclated  noise,  s^t)  is  the  beam. 

A  rough  estimate  of  the  coda  suppression  capabili¬ 
ties  of  the  beam  and  mixed  signal  processor  may  be  had 
by  analyzing  the  results  shown  in  Figure  1  (Shumway, 

1972).  Here  two  signals  recorded  separately  are  super¬ 
imposed  to  simulate  interfering  arrivals.  Signal  2, 
which  in  our  case  represents  the  explosion,  is  shown 
arriving  first.  This  signal  is  actually  a  Tonga  Islands 
earthquake.  Signal  l.the  masking  earthquake,  is  a  Fox 
Islands  earthquake.  Normalizing  all  measurements  to  the 
second  (undisturbed)  cycle  of  the  Tonga  Islands  event,  we 
see  that  the  leakage  of  signal  1  into  the  beam  at  point 
(B)  in  the  figure  is  only  5  db  down  from  the  unprocessed 
signal  shown  at  (A),  while  the  mixed  signal  processor 
has  reduced  signal  1  at  (C)  by  13  db  compared  to  the 
signal  at  (A).  Thus,  in  this  example, leakage  of 


-  4  (K> 


Figure  1.  Least-square  technique  fcr  estimation  of 
presence  of  another  (example).  After  Shumwav  (1972) 


signal  1  into  our  mixed  signal  estimate  is  reduced  by  8  db 
over  the  leakage  observed  in  the  beam  trace. 
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DATA  ANALYSIS 


A  suite  of  five  earthquakes  (Table  1)  recorded  at 
TFO  (Figures  2  and  3,  Table  II)  provides  the  data  base 
for  this  study.  The  filtered  traces  (the  filter  has  0.4 
and  3.0Hz  3  db  points)  are  shown  in  Figures  4  through  8. 
The  7-  and  19-element  subarrays  to  be  examined  consist 
of  elements  Z1  through  Z7,  and  of  Z1  through  Z19, 
respectively.  The  first  arrivals  for  all  the  events 
exhibit  signal- to-noise  ratios  greater  than  20  db. 

Signals  with  high  signal - to-noise  ratios  are  required 
to  insure  that  background  noise  does  not  enter  into 
the  computations. 

For  a  given  event  and  subarray,  the  analysis 
procedure  is  as  follows:  the  onset  of  the  event  (or 
other  common  point  within  the  first  second  of  onset) 
is  picked  on  each  channel.  With  these  common  picks  as 
reference  points,  the  traces  are  shifted  relative  to 
one  another  to  simulate  the  arrival  of  the  event  from 
a  given  distance  and  azimuth. 

Using  the  shifted  traces,  the  beam  and  mixed-signal 
estimates  are  then  computed  for  a  signal  presumed  to 
originate  at  a  selected  test  site.  These  signal  estimates 
represent  the  leakage  from  signal  1  (the  earthquake) 
into  the  estimates  for  signal  2  (the  explosion).  The 
power  for  the  leakage  signals,  when  ratioed  to  the 
averaged  power  for  the  earthquake  arrivals,  yields  a 
set  of  coda  attenuation  figures  for  the  two  processors. 
That  is,  let  R  be  a  measure  of  the  average  rms  input 
signal  power: 
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TABLE  I 

Earthquake  Source  Data 
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Figure  2.  Geographical  location  of  TFO. 
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Figure  3.  Array  configuration  for  TFO 


TABLE  II 

Coordinates  of  TFO  Instruments 
Z1  -  Z19 
SHORT  PERIOD 

LATITUDE  LONGITUDE  ELEVATION  (Meters) 
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2. 2670W 

1491.9 

Z-6 

34 

14 

S5.6970N 

111 

20 

7 .00301V 

1509.4 

1-1 

34 

17 

9.1830N 

111 

21 

24.5490W 

1403.2 

Z  -  8 

34 

21 

42.2900N 

111 

20 

23. 5360W 

1805.4 

Z-9 

34 

21 

10.0020N 

111 

17 

12.9640W 

1569.4 

Z-10 

34 

20 

52.2090N 

111 

14 

4.9130W 

1658.5 

Z-ll 

34 

18 

4. 86  SON 

111 

12 

23. 8630W 

1904.9 

Z- 12  - 

34 

15 

49.5070N 

111 

11 

44. 7380W 

1528.29 

Z  - 13 

34 

13 

48.4020N 

111 

13 

48.3160W 

1513.74 

Z  - 14 

34 

11 

41.4430N 

111 

16 

30.2940W 

1534.0 

Z  - 15 

34 

12 

8.0760N 

111 

19 

8.4440W 

1487.4 

Z-  16 

34 

12 

32.0810N 

111 

22 

13.9550W 

1462.9 

Z  - 1 7 

34 

15 

3.9540N 

111 

23 

28.2740W 

1426.4 

Z-  18 

34 

17 

40.8760N 

111 

24 

39.9550W 

1664.9 

Z- 19 

34 

19 

39.8540N 

111 

22 

32 . 40 10W 

1588.55 

,  N  ,  T  7  1/2 

*  -  h  J  It  I  ^(t)] 

1*1  t*l 


where 

s(t)  is  the  amplitude  of  the  it^1  channel  at  time  t 

T  is  the  length  of  the  sample  analyzed  (the  first 
35  seconds  of  the  signal  in  this  work); 

N  is  the  number  of  channels  processed  (7  or  19  in 
this  report). 

If  §g(t)  is  the  estimate  for  the  leakage  signal  deter¬ 
mined  for  the  beam,  then 

i  T  o  1/2 

-  It  l 

t*l 

is  a  measure  of  the  rms  power  in  §g(t).  Similarly,  a 
measure  of  the  rms  power  in  the  mixed-signal  estimate 
§s(t)  is: 


The  coda  attenuation  figures  are  then  given  by: 
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The  analysis  is  performed  at  a  set  of  uniformly 
spaced  azimuths,  yielding  polar  plots  of  the  coda 
attenuation  capabilities  for  the  beam  and  the  mixed- 
signal  processor.  Note  that  this  analysis  method  is 
equivalent  to  looking  at  the  coda  suppression  obtain¬ 
able  on  that  portion  of  signal  1  which  arrives  ahead 
of  signal  2. 

As  the  events  ranged  from  38  to  85°  distance  from 
TFO,  a  nominal  epicentral  distance  for  the  earthquake 
of  60°  was  chosen  for  this  study.  Further,  computations 
are  performed  at  an  azimuthal  interval  of  30°  (0°,  30°, 
60° , . . . 330° )  ,  with  additional  analyses  performed  every 
10°  in  the  vicinity  of  the  test  site  azimuth.  The  test 
site  used  is  Sen.^palatinsk  (50°N,  78°E),  at  an  epicentral 
distance  of  9S.9°  from  TFO  and  a  station- to- test  site 
azimuth  of  354°.  This  test  site  is  the  epicenter  for 
event  2,  the  (non-existant)  signals  of  which  we  assume 
are  masked  by  the  signals  from  event  1. 
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RESULTS 


Polar  plots  of  coda  attenuation  are  given  in 
Figures  9  through  18.  Results  for  the  7-  and  19-element 
subarrays  are  paired  by  event.  With  but  one  minor 
exception  (Fiji  Islands,  19-element  array,  180°  azimuth), 
the  coda  attenuation  capability  of  the  mixed-signal 
processor  exceeds  that  of  the  beam.  While  up  to  14  db 
improvement  over  the  beam  is  obtained  (Easter  Island, 
19-element  array,  30°  azimuth),  the  improvement  is 
generally  on  the  order  of  3  to  5  db  for  both  subarrays. 
The  maximum  coda  attenuation  obtained  using  the  mixed- 
signal  processor  and  the  19-element  subarray  is  22  db 
(Fox  Islands,  19-element  array,  180°  azimuth),  with  a 
nominal  maximum  value  of  18  db  being  fairly  representa¬ 
tive.  With  respect  to  the  7-element  subarray,  the 
maximum  attenuation  obtained  using  the  mixed-signal 
processor  is  about  17  db  (Fox  Islands,  120  azimuth), 
with  14  db  being  a  more  representative  figure.  In 
general  the  coda  attenuation  obtained  using  the  7-element 
subarray  and  the  mixed-signal  processor  is  comparable 
to  that  obtained  using  the  19-element  array  and  the 
beam  processor. 

To  determine  the  effect  of  narrowing  the  signal 
bandpass,  data  for  the  Fox  and  Tonga  Islands  events 
were  filtered  from  0.8  to  2  Hz  (3  db  points),  and 
analyzed  using  the  beam  and  mixed-signal  processor. 

The  results,  shown  in  Figures  19  through  22,  are 
essentially  the  same  as  those  obtained  previously  using 
a  filter  bandpass  of  0.4  to  3.0  Hz; thus,  there  appears 
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180° 


Figure  9.  Coda  attenuation,  19-element  subarray,  Easter  Island 
event,  seismograms  bandpass  filtered  with  3  db  points  at  0.4 
and  3.0  Hz. 
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Figure  10.  Coda  attenuation,  '/-element  subarray,  Easter  Island 
event,  seismograms  bandpass  filtered  with  3  db  points  at  0.4 
and  3.0  Hz. 


1BD' 


Figure  11.  Coda  attenuation,  19-element  subarray,  Fiji  Islands 
event,  seismograms  bandpass  filtered  with  3  db  points  at  0.4 
and  3.0  Ht. 
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Figure  12.  Coda  attenuation,  7-element  subarray,  Fiji  Islands 
event,  seismograms  bandpass  filtered  with  3  db  points  at  0.4 
and  3.0  Hz. 
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Figure  13.  Coda  attenuation,  19-element  subarray ,  Fox  Islands 
event,  seismograms  bandpass  filtered  with  3  db  points  at  0.4 
and  3.0  Hz. 


Figure  14.  Coda  attenuation,  7-elemeut  subarray,  Fox  Islands 
event,  seismograms  bandpass  filtered  with  3  db  points  at  0.4 
'  ‘  Hz. 
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Figure  IS.  Coda  attenuation,  19-element  subarray,  Hokkaido 
event,  seismograms  bandpass  filtered  with  3  db  points  at  0  4 
and  3.0  Hz.  * 
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Figure  17.  Coda  attenuation,  19-element  subarray,  Tonga  Island 
event,  seismograms  bandpass  filtered  with  3  db  points  at  0.4 
and  3.0  Hz. 


Figure  19.  Coo*  attenuation,  19-element  subarray,  Fox  Islands 
event,  seismograms  bandpass  filtered  with  3  db  points  at  0.4 
and  3.0  Hz. 
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Figure  20.  Coda  attenuation,  7-element  subarray,  Fox  Islands 
event,  seismograms  bandpass  filtered  with  3  db  points  at  0.4 
and  3.0  Hz. 
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Figure  22.  Coda  attenuation,  7-element  subarray, 
Tonga  Islands  event,  seismograms  bandpass  filtered 
with  3db  points  at  0.8  and  2.0Hz. 


to  be  no  advantage  in  using  narrow-bandpassed  data  for 
coda  attenuation. 

The  general  shape  of  the  coda  attenuation  curves 
suggests  that  these  data  can  be  represented  analytically 
by  the  following  formula: 


p  *  a  -  b*cos(6+ot)  db 


(4) 


where  a  is  a 
tion  pattern. 
Least  square 


fixed  rotation  angle  to  orient  the  attenua 
In  the  cases  examined  in  this  work,  a  =  6 
fits  to  the  data  sets  shown  in  Figures  9 


through  18  yield  the  estimates  for  a  and  b  shown  in 
Table  III.  The  average  coda  attenuation  characteristics 


are  shown  in  Figures  23  and  24.  Here,  the  coda  attenua¬ 
tion  curves  have  been  rotated  by  6°,  implying  that  the 
event  to  be  masked  is  arriving  from  the  north. 


While  it  is  convenient  to  describe  the  attenuation 
data  in  the  functional  form  [a  -  b*cos(6+a)  db],  these 
curves  are  only  partially  useful  in  describing  the 
capabilities  of  a  given  array  and  processor.  We  would 
like  to  have  more  general  analytical  expressions  for 
coda  suppression.  Such  curves  can  be  derived  by 
plotting  the  attenuation  data  as  a  function  of  the 
difference  in  the  ray  parameters  for  the  two  signals: 
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Figure  23.  Average 


coda  attenuation,  19-element  subarray 


Figures  25  through  28  display  the  coda  attenuation 
data  as  a  function  of5(dt/dA).  Included  here  are  the 
results  obtained  when  the  source  region  is  taken  as 
the  Tonga  Islands  (A  =  78°)  and  the  masking  event 
(Fox  Islands  event)  is  again  taken  at  60°  distance 
(Figures  29  and  30;  Table  IV).  The  latter  data  set  is 
included  to  help  define  trends  in  the  attenuation 
characteristics,  and  is  not  used  in  subsequent  numerical 
analyses . 

We  are  guided  in  our  choice  of  an  analytical 
expression  for  the  more  generalized  attenuation  curves 
by  the  following  observations: 

I.  Knowing  that  the  generalized  curves  must  pass 
through  (0,  0)  (no  coda  attenuation  possible  when  both 
events  have  the  same  epicenter),  the  data  suggest  that 
for  6  (dt/dA)  small, the  coda  attenuation  curves  are 
linear. 

II.  For  5 (dT/dA)  large,  the  data  appear  to  level 
off  as  if  to  approach  asymptotically  some  limiting 
value. 

Note  that  neither  of  these  conditions  is  satisfied 
by  the  [a  -  b»cos  (9+a)]  model,  which  can  only  describe 
a  portion  of  the  more  generalized  curve. 

One  expression  which  can  be  used  to  describe  the 
attenuation  characteristics  we  observe  is: 

p  .  I  - -  db, 

1.0+5X  +  eX 


(5) 
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Figure  26.  Attenuation  as  a  function  of  the  difference  in  ray 
parameters,  19-element  subarray,  mixed-signal  processor. 
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Figure  27.  Attenuation  as  a  function  of  the  difference  in  ray 
parameters,  7-element  subarray,  beam. 


Figure  28.  Attenuation  as  a  function  of  the  difference  in  ray 
parameters.  7-element  subarray,  mixed-signal  processor. 
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Figure  30.  Coda  attenuation,  7-element  subarray,  Fox  Islands 
event,  seismograms  bandpass  filtered  with  3  db  points  at  0.4 
and  3.0  Hz.  (Signal  2  from  Tonga  Islands). 


TABLE  IV 

Coda  Attenuation  Parameters 
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where 


X  =  6(dt/dA) 

Using  equation  (5)  we  obtain  the  least  square 
solutions  shown  in  Table  V.  These  solutions  are  shown 
graphically  in  Figures  25  through  28.  While  some  of 
the  analytical  expressions  do  not  behave  well  at  very 
large  values  of  X,  all  of  the  solutions  are  acceptable 
in  the  region  of  interest.  Note  that  for  either 
array,  the  value  for  3  determined  from  the  mixed-signal 
data  is  almost  six  times  greater  than  that  determined  for 
the  beam  data.  As  3  is  the  slope  of  the  attenuation 
curve  for  X  small,  this  result  demonstrates  the  high- 
reoolution  capability  of  the  mixed  signal  processor 
over  that  of  the  beam. 

Solutions  to  equation  (5)  provide  considerably 
better  fits  to  the  data  than  do  those  solutions  based 
on  equation  (4).  We  therefore  recommended  employing  the 
more  generalized  solutions  in  determining  the  number  of 
opportunities  for  masking  the  signals  from  an  explosion 
in  the  coda  of  an  earthquake.  When  computing  the 
number  of  opportunities,  a  given  number  of  stations  in 
a  worldwide  monitoring  network  would  be  given  coda 
suppression  capabilities  based  on  the  number  of  sensor 
elements  specified  (7  or  19)  and  the  type  of  array 
processing  to  be  performed  (beam  or  mixed-signal 
processor).  That  is,  a  given  station  would  be  assigned 
one  of  the  coda  suppression  curves  listed  in  Table  V. 
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TABLE  V 

Coda  Attenuation  Parameters 


where  X  =  & 


19 -ELEMENT 

SUB ARRAY 

BEAM 

mixed-signal 

PROCESSOR 

0.98 

S.74 

0.14 

0.91 

0 

0.31 

0.006 

0.43 

3X  +  yX2 
+  6X  +  eXZ 


7 -ELEMENT  SUBARRAY 


BEAM 

mixed-signal 

processor 

0.37 

2.21 

0.04 

0.69 

0.00 

0.53 

0.00 

0.005 

^I'Lb  - 


An  attenuation  factor  based  on  the  vector  difference 
in  the  ray  parameters  for  the  two  epicenters  under 
investigation  can  then  be  determined,  and  the  coda  of 
the  masking  earthquake  scaled  accordingly. 

As  a  final  note,  let  us  compare  the  coda  attenua¬ 
tion  obtained  for  a  wide-aperture  7-element  TFO 
subarray  (Zl,  Z9,  Zll,  Z13,  Z15,  Z17,  and  Z19)  with 
that  obtained  previously  for  the  19-elem'nt  and  small 
aperture  7-element  subarrays  (Table  VI,  Figures  31  and 
32)  The  data  analyzed  are  from  the  Easter  Island  event. 
For  the  wide-aperture  7-element  subarray  and  with  the 
exception  of  those  azimuths  for  which  the  difference  in 
the  ray  parameters  is  small,  the  beam  and  mixed-signal 
processor  perform  about  equally  well;  the  mixed-signal 
processor  yields  only  1  to  2  db  additional  attenuation 
over  the  beam.  Further,  the  coda  attenuation  obtained 
using  the  wide-aperture  7-element  subarray  and  the 
mixed-signal  processor  is  about  the  same  as  that 
obtained  for  all  19-elements  and  the  beam,  as  well  as 
that  obtained  for  the  7  inner  elements  (Z1-Z7)  and  the 
mixed-signal  processor.  Finally,  for  small  ray-parameter 
separation,  the  wide-aperture  7-element  subarray  and 
the  mixed-signal  processor  produce  about  the  same 
attenuation  as  does  the  19-element  subarray  and  the 
mixed-signal  processor  (-9  db). 

In  sum,  a  wide  aperture  array  and  the  mixed-signal 
processor  are  more  important  than  the  number  of  elements 
for  significant  coda  attenuation  when  the  difference  in 
the  ray  parameters  is  small.  When  the  ray  parameter 
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Figure 
event , 


31.  Comparison  of  coda 
19-element  subarray  and 


attenuation  figures,  Easter  Island 
wide- aperture  7-element  arrav. 
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Figure  32.  Comparison  of  coda  attenuation  figures,  Easter  Island 
event,  small-aperture  7-element  subarray  and  wide-aperture 
7-element  subarray. 
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difference  is  large,  and  when  using  the  mixed-signal 
processor,  the  wide-aperture  7-element  subarray  performs 
about  as  well  as  the  small-aperture  7-element  subarray. 
Where  using  the  beam,  however,  the  wide-aperture  7-element 
subarray  performs  substantially  better  than  does  the 
small-aperture  7-element  subarray,  the  improvement  being 
on  the  order  of  5  db. 
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CONCLUSIONS 


A  comparison  has  been  made  of  the  coda  attenuation 
capabilities  of  the  beam  and  mixed-signal  processor 
using  7-  and  19-element  TFO  subarrays.  The  major  source 
region  was  taken  as  Semipalatinsk  (96°  distance  from 
TFO),  while  the  interfering  signals  from  various 
earthquakes  were  time-shifted  to  simulate  arrivals  from 
various  azimuths  at  60°  distance.  Our  major  conclusions 

are : 

1.  The  coda  attenuation  capability  of  the  mixed- 
signal  processor  exceeds  that  of  the  beam.  While  up  to 
14  db  improvement  was  obtained,  the  improvement  is 
generally  on  the  order  of  3  to  5  db  for  both  subarrays. 

2.  A  nominal  value  of  18  db  is  representative  of 
the  maximum  coda  attenuation  obtained  using  the  mixed- 
signal  processor  and  the  19-element  array.  For  the 
7-element  array,  the  corresponding  figure  is  14  db. 

3.  The  coda  attenuation  obtained  using  the  7-element 
sub array  and  the  mixed-signal  processor  is  comparable 

to  that  obtained  using  the  19-element  array  and  the 

beam. 

4.  No  advantage  was  observed  in  using  data 
filtered  in  the  band  0.8  to  2.0  Hz  over  that  data  which 
was  filtered  in  the  band  0.4  to  3.0  Hz. 

5.  Representative  TFO  coda  attenuation  curves  for 
the  mixed-signal  processor  are: 
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19  elements: 


p  =  5.74X  ♦  0.91X2  db 

1  +  0.31X  +  0.043X2 

7  elements: 

p  =  2.21X  ♦  0 . 69X2  db 

P  l  +  0.53X  +  0.005X2 

where 

x  5 

6.  For  X  small,  the  mixed-signal  processor  is 
six  times  more  sensitive  to  changes  in  X  than  is  the 
beam. 

7.  Preliminary  investigations  at  TFO  using  a  wide- 
aperture  7-element  array  suggest  that  for  small 
differences  in  the  ray-parameter  vectors,  a  large- 
aperture  array  and  the  mixed-signal  processor  are  more 
important  than  the  number  of  elements  for  significant 
coda  attenuation  (~9  db). 
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